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Abstract 

We propose a scheme for nonlinearizing linear equations to generate integrable nonlinear systems of 
both the AKNS and the KN classes, based on the simple idea of dimensional analysis and detecting the 
building blocks of the Lax pair. Along with the well known equations we discover a novel integrable 
hierarchy of higher order nonholonomic deformations for the AKNS family, e.g. for the KdV, the 
mKdV, the NLS and the SG equation, showing thus a two-fold universality of the recently found 
deformation for the KdV equation 6kdv. 

I. INTRODUCTION 

In spite of the significant achievements in the theory of integrable systems over the last hfty years, 
it remains a mystery: what is integrable nonlinearity, i.e. how to identify a priory those nonlinear terms 
which make an equation integrable. 

Nevertheless one notices that, although integrable equations can be of diverse types, the integrability 
filters out the nonlinearity to some specific forms, as evident from the well known equations like the 
nonlinear Schrodinger (NLS) equation, the derivative NLS (DNLS), the Korteweg-de Vries (KdV) equa- 
tion, the modified KdV (mKdV), the sine-Gordon (SG) equation etc.p]. Moreover, all such integrable 
nonlinear PDEs exhibit exact soliton solutions, localized and stable form of which is believed to be due to 
a fine balance between their linear dispersive and nonlinear terms. Therefore it may not be unexpected 
to expect that a linear PDE, the dispersive part of which would balance the nonlinear term, should have 
the information of the integrable nonlinearity hidden in it. Consequently, there must be a scheme for 
nonlinearizing linear equations for generating integrable systems. However attacking this anti-intuitive 



question directly seems to have been avoided in the hterature, although the Ablowitz-Kaup-Newell-Segur 
(AKNS) scheme [1] and the celebrated Sato construction [2] aimed around this line and a recent attempt 
came close to it [3|- 

We propose here a direct but easy scheme for nonlinearizing the linear equations 

iQt = Qxx, Qt = Qxxx, (^xt = 0, (1) 

etc. together with their conjugates, which are in fact the linearized parts of the well known equations like 
the NLS, the DNLS, the KdV, the SG etc. Our aim is to generate first the known integrable systems, 
proving thus the effectiveness of our alternative scheme based on a simple physical idea of dimensional 
analysis and then, continuing the process search for new integrable nonholonomic deformations for all 
known equations. 

It is evident that the direct nonlinear extension of a linear equation is not unique, since it alone can 
not identify the integrable nonlinearity and one needs some extra structure for filtering out the integrable 
cases. Our strategy therefore is to stick from the beginning to a Lax pair, a property inherent to all 
integrable systems, by defining a pair of naive Lax operators for the linear system itself and 

then nonlinearize them to a genuine pair ([/, F), which finally generate the nonlinear integrable equation 
through their flatness condition 

Ut-V^ + [U, V] = (2) 



It is true that several sophisticated formal methods like the prolongation technique, the Painleve trun- 
cation [1], the Sato theory [2j, the AKNS method p!] etc. are available in the literature for constructing 
the Lax pair of a nonlinear integrable equation. However our alternative construction, based on a fun- 
damental concept of dimensional analysis, is perhaps the simplest and the most physical one. Exploiting 
only the scaling dimensions and identifying the constituent building blocks of the Lax operators, hidden 
already in the linear equations, we can construct a Lax pair quite easily to any higher order for both 
the AKNS p] and the Kaup-Newell (KN) |;5j spectral problems. Interestingly, by regulating the scaling 
dimension of the field we can generate uniquely both the NLS and the DNLS equations, starting from the 
same linear Schrodinger (LS) equation. Finally as an important application we discover new integrable 
hierarchies of nonholonomic deformations for the well known equations like the KdV, the mKdV, the 
NLS, the SG etc. The nonholonomic deformation in this field theoretical context is given by some differ- 
ential constraints on the deforming function. A recently discovered 6-th order integrable KdV (6KdV) 
equation, attracting much attention [El [71 [8], happens to be the lowest order deformation of this kind 
for the KdV equation. We are able to generalize here not only the earlier result on the deformed KdV 
and the 6KdV by unveiling a novel hierarchy of their higher order deformations, but also to show the 
universality of such perturbations, by discovering integrable nonholonomic deformations for all members 
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of the AKNS family, including their higher order integrable hierarchies, a task termed as highly desirable 
but quite impossible in a recent study [YJ. Interestingly, such deformed integrable systems can be rep- 
resented also as known nonlinear equations with additional perturbative terms subjected to differential 
constraints. Such perturbations, contrary to the usual belief, not only preserve the integrability of the 
original equation but also yield richer properties. 

Nonlinear integrable systems like the NLS, the mKdV and the SG equations with various pertur- 
bations appear in many physical situations, e.g. fluxons in Josephson junction, parametrically driven 
damped molecular chain, nonlinear optical fiber communication, ferromagnet in variable magnetic field, 
nonlinear Faraday resonance etc. [9j. However in almost all these perturbed systems the integrability - 
the most cherishable property of a nonlinear system - is usually lost. As a result all their related pre- 
cious properties like higher conserved quantities, exact soliton solutions, elastic scattering of solitons etc., 
enormously useful in physical applications, are also lost. Therefore constructing integrable perturbed 
systems, preserving their complete integrability and exact soliton solutions, is a challenging physical 
problem, which we are able to solve here to certain extend. 

Our integrable perturbed equations exhibit the usual exact N-soliton solutions, but with an unusual 
accelerated or decelerated motion. Recall that the standard solitons move with a constant velocity and 
behave like particles undergoing elastic collisions with other solitons, which have been of significant 
importance and application in various physical problems. In the present context of integrable perturbed 
equations the accelerated solitons behave like particles under external forces. They are subjected to 
perturbations, which are themselves solitonic in nature. In fact the time-dependent asymptotic value of 
such a perturbing function acts like a force sitting at the space boundaries and controlling the motion of 
the field soliton, depending on the nature of which the soliton can accelerate or decelerate. 

Accelerated solitons appear in many physical models like inhomogeneous plasma [lOj , information 
transfer in the DNA chain [llj , transport of solitons and bubbles [13] etc. Therefore the present integrable 
models with exact accelerating solitons showing elastic soliton scattering could find applications in 
practical situations. 

The arrangement of the paper is as follows. Sec. II introduces our nonlinearization scheme. Gen- 
eration of specific nonlinear integrable systems is presented in Sec. III. New integrable nonholonomic 
deformations of the AKNS system are detailed in Sec. IV, with the concrete examples given in Sec. V. 
In Sec. VI the present result is listed and compared with the existing ones. Sec. VII is the concluding 
section followed by the bibliography. 

II. THE NONLINEARIZATION SCHEME 

Combining field q and its conjugate r in a matrix as U^^^ = qa^ + ra^ , where = ± ^o"^) 
and a"" ,a = 1, 2, 3 are standard 2 x 2-Pauli matrices, we can clearly express our starting linear dispersive 
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equations ([T|) as iU^ = asUxx , = UxJx-, etc. These linear equations in turn can be expressed by 

easy observation as a linear flatness condition U^^^ — V^^^ = 0, which is defined by ignoring the nonhnear 
commutator term between Lax operators. This naturahy identifies the pairs of naive Lax operators 
V^^^), i = 2,3, . . . for the consecutive higher order hnear equations as 

^7(0 = iXas + iU^'\ and V^'^ = a^U^'l V^^ = iU^^J , (3) 

etc. with trXjW = trVf'' = 0, i = 2,3, since trcja = trU(°) = tr(o-3U(°)) = 0. The significance of 
parameter A, which enters here only as an additional constant and can be ignored without any loss, will 
become clear in the process. 

The Next important point is to detect the scaling dimension of the Lax operators as well as to identify 
their crucial building blocks (BB), existence of which seem to have been overlooked in most of the work 
including the AKNS construction. 

A. Scaling dimension and building blocks of the Lax pair 

We consider length L as the fundamental dimension and from the linear evolution equation with 
A^-th order dispersion: qt = qxx . . . x read out the dimension [•] of time as [T] = . The Lax operators, 

JV 

generating infinitesimal space and time translations, as evident from the Lax equations = U^, = 
V^, naturally should have the scaling dimensions (SD) as [U] = L^^, \Vn] = [T]~^ = . Assuming 
that the linear Lax pair U^^\V^^ share the same scaling dimensions as well as the same building blocks as 
their nonlinear counterpart, we find from the explicit structure ([3]), the SD of their constituent elements 
as [A] = \U^^'>] = and at the same time identify the building blocks (BB) for both the Lax operators 
as the set { A, f/^*^) } and its x-derivatives: U^\u^x ■, etc. , which enter in their construction linearly. 
Our crucial conjecture is that these linear Lax pair can be nonlinearized by adding all possible nonlinear 
combinations of the same BB, guided by the dimensional argument. Note that, the condition \U^^^] = L^^ 
fixes the SD for the field in this case as [q] = [r] = L^^. However, since a linear equation as such can not 
put any restriction on the SD of its fields, this input is not unique and we show below that a different 
fixation of the SD for the field would lead to another construction of the Lax pair. 

Recalling that a nonlinear variable change can generate from a given integrable system other gauge 
equivalent equations [H], we concentrate here only on fundamental integrable equations like NLS, KdV 
DNLS etc. belonging to the AKNS or the KN family, while the gauge equivalent systems can be obtained 
trough simple transformations. 

III. GENERATION OF NONLINEAR INTEGRABLE SYSTEMS 

The dimension of the field [q] and its conjugate [r] as well as [A], which are not fixed by the starting 
linear equations ([1]), should be given in our nonlinearization scheme as input. We show that for [q] = 
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[r] = [X] = L ^, we get the AKNS systems, while [q] = [r] = [X] = L 2 yields the KN hierarchy. 

A. AKNS integrable hierarchy 

We consider first the case [q] = [r] = [X] = L~^, which gives through the above construction: [U^^^] = 
[[/(Oj = L^^. We intend to construct now the space-Lax operator U{X) G s/(2),trU(A) = 0, naturally 
with [U{X)] = L-i, out of f/(') by using the BB: {X,U^^^}. Therefore the only possibility left from the 
dimensional argument, which only allows summation of the terms having the same dimension, is to take 
U{X) = = i(Ao"3 + U^^^), recovering thus the well known AKNS Lax operator [1] uniquely. However 
for constructing the corresponding time-Lax operator l^(A) G sl{2) with trV2(A) = and SD: L^'^, in 
addition to V^*''^ = (73Ux^\ nonlinear terms l^*-"'^ with the same SD: L^^ are to be constructed out of 
the same BB {A, U^'^^}, through their nonlinear products and powers like A2, XU^'^\ (C/(°))2, following the 
dimensional argument. Note that since [dx] = , the derivative term U^"^ can appear only in the linear 
part, again from the dimensional analysis. 

Therefore this set of nonlinear terms is sufficient to construct the nonlinear part of the time-Lax 
operator: ¥2^^^ = i{k2X^cj3 + kiXU'^^^ + kQa3{U'^^'>)'^), upto the integer coefficients /c„,n = 0,1,2. Note 
that fixing the dimensionless integers with different terms in V2(A) = ¥2"^ + V^*-""'^ goes beyond the 
scope of the dimensional argument, which however is achieved from the flatness condition of ?7(A), V2(A), 
yielding several consistency relations at different powers of A. Interestingly, the number of such relations 
obtained is just sufficient to determine the numerical coefficients of all the terms without any ambiguity. 
For example, it is easy to check that, the condition ^ for C/(A), V2{X) yields three consistency conditions 
at three different powers of A"', n = 0, 1, 2, of which the condition at n = fixes k^ = —1, that at n = 1 
fixes the integer ki = 2, while n = 2 gives the relation k2 = ki. Thus all the integer coefficients are 
obtained exactly, determining the structure of the nonlinear part unambiguously as 

= 2iAV3 + 2U[/W - icT3(C/(°))^ (4) 

with trV2°^^ = 0, which constructs the complete time-Lax operator uniquely as V2(A) = V^*''^ + V^^"'^ . 
Note that the asymptotic solutions of the Lax equation giving e^*^^ clarifies the physical meaning of 
parameter A as the momentum and justifies the choice of its SD: L^^. At the final step we obtain the 
integrable nonlinear equation we are searching for, again from the flatness condition, where a nonlinear 
term 2as{U^^^)^ appears now in addition to the initial linear equation iuj^^^ = a^Ux^ . For r = ±q* 
this nonlinear term reduces to ib2|gpg, which together with the starting LS equation yields finally the 
integrable NLS equation, completing our nonlinearization process. 

The next higher order linear equation qt = qxxx follows a similar procedure in the nonlinearization 
process. The space-Lax operator being the same, one has to build only the time-Lax operator V3(A) by 
adding to the naive operator V^^'^ = iU^x a nonlinear part VJ^'^''' with SD L~^, which is to be constructed 
again from the same BB : {A, U^^^} and its derivative following our conjecture. Observe that the 
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derivative term can appear now in the nonlinear part on dimensional ground. By all possible nonlinear 
combinations as powers and products of the BB with total SD = L^^ and maintaining trV3(A) = 0, we 
can, similar to the above, construct uniquely : 

yielding finally V3(A) = V3'^ + V^^^\ The flatness condition of the pair U{\), V3(A) fixes again the integer 
coefficients in different terms of ([5]) and yields the integrable equation, by adding only one nonlinear term 
([/(o)^2^_(o) ^YiQ starting linear dispersive equation ul^^ = U^Jx ■ For r = k = const., q = u the 
nonlinear term reduces to k6uUx, while for r = q = v to Qv'^Vx, generating thus the well known integrable 
KdV and mKdV equations, from the third order linear equation we started with. 

Similarly one can continue building the hierarchy of integrable equations, starting from the arbitrary 
higher order linear equation qt = qxx...x and nonlinearize it following the procedure as above and using 
the same BB and similar dimensional argument, as we have conjectured. Since the space-Lax operator 
U{X) remains the same, the task is to construct only the time Lax operator Vn{X) with the SD for 
arbitrary N, out of the same BB: {A, U^^^} by all possible combinations like 

{x)\^)\u^o^r (6) 

with the dimensional constraint k + l + m = N. It is easy to check that this partitioning of N determines 
the total number of terms that can appear in Vat (A) as L = 1 + ^N{N + 1), which interestingly tallis 
with our result obtained above, giving L = 4 for the NLS with N = 2 and L = 7 for the KdV and mKdV 
with = 3. Thus we can estimate the exact number of terms that should appear in any higher order 
Lax operator, without even calculating their explicit form. This unique feature of our nonlinearization 
scheme, obtained from the dimensional argument and the identification of the BB, is absent in all other 
available methods. 

B. Comparison with the AKNS scheme 

It is true that, though we have presented a simple and original scheme for nonlinearizing linear 
equations to integrable systems, based on the dimensional analysis and using the building blocks of the 
Lax pair, we could reproduce so far only the known integrable equations, which are obtainable also 
through the AKNS scheme In-spite of this fact and before presenting our new result, we intend to 
show that in comparison with the AKNS proposal [1], the construction of time-Lax operator y/v(A), vital 
in generating integrable equations, is significantly simpler in the present scheme. 

Recall firstly, that in the AKNS method one starts with a given nonlinear field equation and not from 
its linearized version, as done here. 

Secondly, for obtaining the explicit form of Vm{^) in the AKNS method, one has to expand it in 
the powers of spectral parameter: Vn{X) = J2n=i — X^V^'^\ where y("),n = 1,2, .. . are number of 
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2x2 unknown matrices with 4:N number of unknown functions, some of them being complex. In our 
construction on the other hand, the unknown coefficients are only N in number and moreover they are 
only real integers. This is because we have already identified the building blocks of the Lax operators 
made up from parameter A, a known matrix U^^^ and its derivatives, and the dimensional argument, 
which is our key ingredient, guides us to collect them in the right combination. 

Thirdly, in the AKNS scheme the 4A^ unknown functions are to be determined by solving different 
partial differential equations. Our N unknowns being only real integers are obtained from simple algebraic 
equalities. 

This simplicity of our scheme becomes more evident for higher values of N, where Vat contains 
L = 1 + ^N{N + 1) number of terms. Note that determining a priori the number of terms that should 
be present in any Vm, as we have shown above is again easy in our scheme but not in the AKNS. 

Fourthly, while for the KN spectral problem the AKNS type direct expansion becomes even more 
complicated, our scheme covers the KN system in an unified way. Just a change in the scaling dimension 
of the field switches our scheme from the AKNS to the KN family, as we demonstrate below. 

Finally, as an important application we discover new integrable hierarchies of nonholonomic deforma- 
tions for all members of the AKNS family (see Sec. V), systematic account of which is absent in AKNS 

C. KN integrable hierarchy 

Our nonlinearization scheme , as we show here, is universally applicable for generating the KN hier- 
archy, which is not readily available in the literature due to its apparent complicacy. 

We follow again the same line of argument, making a slight change in the input of the SD of the field 

[q] = [r] = [A] = L^2. This small deviation however dramatically changes the outcome and resolves 
the puzzle we faced in generating uniquely two different nonlinear equations, namely the NLS and the 
DNLS equations, starting from the same linear LS equation. 

Note that the SD of the linear Lax operator has been changed now to = L"2, while that for 

the genuine space Lax operator must always be [f7(A)] = L^^. Therefore the nonlinearization induced 
by the dimensional argument should give U{X) = AC/^'^ = iX'^cr^ + iXU^^\ which reproduces correctly the 
well known KN Lax operator [5j. Note that, in principle, one can also add on dimensional ground, a 
term like icr^iU^^'^)'^ to this Lax operator. However, such an addition does not produce any independent 
equation that are not derivable from those obtained without its addition. 

For constructing the corresponding time-Lax operator V]\[{X) through our nonlinearization, one should 
remember that the situation is a bit different here, since though the dimension of Vn{X) must remain as 
L"^, the SD of its BB: {A, C/(°)} has been changed. 

Let us consider first the case with N = 2 and start again from the linear Schrodinger equation. 
We repeat the above procedure for the NLS case, but with the present change in the dimension of the 
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BB, which should lead to the construction V2{X) = AV^ + ^''"''^(A) due to the dimensional constraint 
where V^^^ = a^U^^^ is the naive linear operator as in ([3]). The nonlinear part y^"'^(A) should be 
constructed therefore by nonlinear combinations of the BB with dimension 

L-2 = [A][(C/(°))3] = [A2][(?7(°))2] = [X\U^^'^] = [X\ (7) 

Note again that though terms like zcj3({7^'^^)'^, \U^^\ U^"^] match in dimensionality and hence can also be 
added, they only yield equations that are derivable from those, obtained without their addition. This 
identification together with V2(A) G s/(2), trV2(A) = 0, constructs the time-Lax operator 

V2{\) = {a^U^^^ + f(C/(°))=^)A - ia^iU^'^^fX^ + 2iC/WA3 + 2ia3A^ (8) 

uniquely , with the integer coefficients fixed as above from the consistency condition. This condition 
generates also the nonlinear equation with the integrable nonlinearity i[{U^^^)'^)xi which for r = itq* 
adds a term ±.2(\q\'^q)x to our starting LS equation iqt — qxx = 0, yielding the DNLS equation, as 
required. 

Similarly the higher order linear equations we considered above nonlinearize now to a different set of 
integrable equations belonging to the KN hierarchy, due to the changed SD of the BB. Without giving the 
details we just mention, that the construction, though a bit tedious, is quite similar and straightforward 
and given again in the form ([6|), with the total number of terms appearing in Vat (A) as L = 1 + N"^ , with 
arbitrary . Note that this number of terms L, resulting from the new SD constraint: \{k + m) + l = N, 
is much higher than the number L appearing in its AKNS counterpart, showing why the KN systems are 
more complicated than those of the AKNS. Check that N = 2 gives L = 5 as we have obtained above for 
the DNLS case. 

IV. NEW INTEGRABLE HIERARCHIES FOR THE AKNS FAMILY WITH NONHOLO- 
NOMIC DEFORMATIONS 

Note that the emphasis in the above nonlinearization scheme is to construct the time-Lax operator 
Viv(A) = V^'"* + V^"'^(A), by building up the nonlinear part involving positive powers of A and confining 
to the elements only from the identified building block U^^^ = qa^ + ra^ , which means to construct 
the nonlinear combinations of the same basic fields q, r and their derivatives. Now we intend to extend 
this construction further by going beyond the conjecture of the BB and introducing new perturbative 
functions with SD L > N. Thus keeping U{X) same as above, we add a nonlinear deformation V^'^'^'^'' (A) to 
the V]sf{X) by including negative powers of A^*^, n = 1, 2, ... in the dimensional argument, which we have 
ignored so far. This simple extension, as we show, would discover a completely new class of integrable 
perturbed equations with hierarchy of higher order nonholonomic deformations. Thus it creates a novel 
two-fold integrable hierarchy [12j by adding a deformation hierarchy to each of the known AKNS hierarchy. 
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The well known equations like the NLS, the KdV, the mKdV and the SG can be deformed by perturbing 
functions with higher and higher order nonholonomic constraints, preserving the original integr ability. 
The present result thus shows a two-fold universality for the nonholonomic deformations, found recently 
for the KdV equation O [71 IS] , since one can cover now the entire AKNS family , originating from this 
single model, and at the same time can discover a new integrable deformation hierarchy for each members 
of this family. 

It should be noted in this context, that the use of negative powers of the spectral parameter in the 
time-evolution operator was considered also in some earlier occasions in the long history of integrable 
systems. However this was either limited [H [15], partial [16j or camouflaged [17j. The novelty of our 
result lies in the fact, that within this extremely well studied field we could discover in a simple way 
a class of new integrable systems with a novel two-fold integrable hierarchy, which can be interpreted 
as perturbed equations with integrable nonholonomic constraints. In this construction of the time-Lax 
operator, we do not confine to the BB, which involve only the basic fields, but include a series of new 
perturbing functions, each influencing the basic field. This situation can simulate an interesting device 
for controlling the basic solitons through multiple intervention, though remaining within the scenario 
of the exact solvability. Such an idea was partially realized in the fiber optics communication through 
doped media [E]. 

For presenting our new equations we start with V2(A) for the AKNS system constructed above and 
deform it first by y('^^^)(A) = ^X'^G^^^ with a matrix function G^^) with trG^^^ = and SD L^^. Since 
our intention is to introduce new perturbing function g into the system, we derive from the flatness 
condition its structure as an integrable deformation of the second-order AKNS equation: 

iqt - qxx - '^{qr)q = g, (9) 

with the perturbing function, given by the matrix element g = ^[2^ subjected to the nonholonomic 
differential constraint: 

G« =i[C/(o),G«]. (10) 

Remarkably, we can include further perturbation into the system, which would interact with the basic field 
as well as with the initial perturbation. This process can also be viewed as putting higher order differential 
constraint on the original perturbation. To achieve this higher perturbation we extend V^'^''J^\X) with 
another deforming term ^\~'^G^'^\ where G^^^ is a matrix function of SD L~^. Integrability condition 
([2|) now leads to a further deformation of Q given by the higher order constraint 

=i[C/W,G«]+i[a3,G(2)], 42) ^. [^(0)^^(2)]^ 

Generating such higher order integrable perturbations can be continued recursively by adding in V^'^^^^ (A) 
more and more terms as ^X~^ G^^\ j = 1, 2, ... n, with arbitrary n. New matrix functions G^^^ have the 
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scaling dimension j + 2 anf trG^j^ = 0. This would result to a new integrable hierarchy of nonholonomic 
deformations for ([9]), given recursively as 

G^^) = i[U^o),G^% (12) 

which clearly reduces to (llOp for n = 1 and to (llip for n = 2. 

Exactly in a similar way we can build up the hierarchy of integrable perturbations for each member 
in the known AKNS hierarchy of higher nonlinear equations, which would lead thus to a novel two-fold 
integrable hierarchy. In one the same nonlinear equation is perturbed by a function with increasingly 
higher order differential constraints and in the other all different higher nonlinear equations are deformed 
by the same perturbing function. 

The most important fact about the hierarchies of all perturbed equations thus generated is that, 
they are completely integrable systems and are exactly solvable by the inverse scattering method (ISM), 
yielding A^- soliton solutions. This follows from the fact that, for constructing such nonholonomic defor- 
mations, we have started from the Lax pair, keeping the space-Lax operator U{X) same as the original 
one, while deforming the time-Lax operator V{X). Therefore the scattering problem, which is central 
to the ISM remains the same, while the time evolution of the spectral data only gets changed for the 
deformed models. Consequently, in such perturbed systems along with the exact soliton solution for the 
basic field we can find the exact solution for the perturbing function, which intriguingly takes also the 
solitonic form. Moreover, one can analytically study the soliton dynamics as well as the scattering of 
multiple solitons, a task impossible to carry out under usual perturbation with a known function. It is of 
significant practical importance that, though in such perturbed integrable systems, the exact ISM is ap- 
plicable for extracting the soliton solution, we can bypass this involved and lengthy procedure and achieve 
the same result by taking the well-known soliton solutions of the undeformed system and deforming them 
by suitably choosing their time-evolution, given by the deformed soliton velocity kunjpa082. 

Remarkably, a perturbed soliton in the present set up behaves like a particle driven by a force and 
exhibits an accelerated or decelerated motion, depending on the nature of the deformation. Note that 
unlike known situations [19] the variable soliton velocity occurs here without any apparent inhomogeneity 
and within the framework of an isospectral fiow. 

V. NEW INTEGRABLE EQUATIONS: SIMPLE EXAMPLES 

Using the above construction suitable for equations belonging to the AKNS spectral problem we can 
now analyze in detail the new class of integrable perturbations for all members of this family, namely the 
KdV, the mKdV, the NLS and the SG equations. In particular constructing the matrix Lax pair for each 
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of them we can find the exphcit form for ah those deformed equations, explore their higher deformations 
and most importantly, obtain their exact N-soliton solutions. 

Along with these perturbed equations we can also study their two-fold integrable hierarchies, as 
developed above. However we present here only the simplest equations in this hierarchy, which are the 
most important cases involving the perturbation of the well known equations and represent the lowest 
order nonholonomic deformation, obtained from (jl2p with n = 1 or 2. 

A. Integrable perturbation of the KdV equation 

Recent findings of the integrable nonholonomic deformation of the KdV equation, equivalent to a 6th- 
order KdV equation [6], which apparently contradicts the accepted notion of nonexistence of any even- 
order equation in the KdV hierarchy, arose considerable interest [318]. However, the exact integrability 
of this system through AKNS type matrix Lax operator or its exact N-soliton solutions through the ISM 
could not be established. Consequently, the integrable hierarchies allowed by this system as well as the 
dynamics of the solitons, including their nature of scattering could not be explored. 

We on the other hand have developed the formalism in Sec. IV for constructing AKNS type Lax 
operators for all of its members with nonholonomic deformations, which should solve completely the 
remaining unsolved aspects of the deformed KdV equation, as we show below. 

Recall that in the AKNS hierarchy KdV type equations appear at = 3 level, for which we have 
derived explicitly the time-Lax operator V3(A), through our nonlinearization scheme. Now we intend to 
deform this Lax operator by an additional ^('^^■^^(A) matrix, which would result to a deformed 3rd-order 
AKNS equation 

qt - Qxxx - Qiqr)qx = 9x, (13) 

with the multiple nonholonomic constraints as in ()12p on the deforming matrices G^^\ having higher 
scaling dimension: j + 3. For deriving the perturbed KdV equation we have to make the standard 
reduction r = 1 , q = u in ([13]) to yield 

Ut - Uxxx - Quua: = gx{t, x). (14) 

For fixing now the lowest order constraint on the deforming function, we have to add two terms with 
n = 1, 2 to V^'^'^-f\X) = ^{X~^G^^^ + A^^G^^-*). It is important to observe that adding only one deforming 
term in this case gives trivial result, while the double-deformation yields from pTj) the explicit structure 

G^^) = -i^gxa^ + {g + c)a~ + iea~^ , = iug^. (15) 

and the same compatibility condition leads to the nonholonomic constraint on the perturbing function 
as 

gxxx + '^ugx + 2ux{g + c) = Q. (16) 
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Therefore the coupled system (jl4til6p gives finally the integrable perturbation of the KdV equation, 
recovering the recent result [6] . We can find now the exact N-soliton solution for this system of equations 
(|16p through the application of the ISM. Referring to [12] for details we present here only the 1- soliton 
solutions for both the KdV field u{x,t) and the perturbing function g{x,t) as 

u{x,t) = -^sech^^, ^ = k{x + vt) + (p, (17) 
gix,t) = c(t)(l - sech^e), (18) 

where v = vq + Vd, with vq being the usual constant velocity of the KdV soliton, while Vd = 
with c{t) = J dtc{t) is its unusual time-dependent part, induced by the deformation. The perturbing 
function, itself taking the solitonic form as (jlSp . drives the field soliton (jl7p through its asymptotic value 
lim|3..|^oQ (7 = c{t). Therefore, in such perturbed integrable systems one can control the motion of the 
soliton from the space-boundaries by tuning an arbitrary time-dependent function c{t), a feature seems 
to be of immense physical significance for practical applications. Fig. la shows the dynamics of this 
perturbed KdV soliton moving with a constant deceleration, obtained for linear function c{t) = cqI with 
Co < 0. The exact 2-soliton solution for this perturbed KdV equation, which can also be derived explicitly, 
shows beautiful elastic scattering of the accelerated solitons (see Fig. lb). 




Figure 1: Dynamics of the exact soliton solution of the integrable perturbation of the KdV equation for 
the field u{x,t). a) 1-soliton with the usual localized form but with an unusual deceleration, refiected 
in the bending of soliton trajectory in the (x,t)-plane, b) 2-soliton with the usual elastic scattering and 
phase shift, but with the dynamics dominated by an unusual accelerated motion. 

B. Integrable perturbation of the mKdV equation 

Integrable equation with a nonholonomic deformation is driven by a perturbing function, which in 
turn is subjected to a differential constraint. Therefore, analogous to the deformed KdV, constructed 
above, we can derive a new deformed mKdV equation again from (jl3p . but with reduction q = r = v. 
Interestingly, in this case a single deformation with n = 1, corresponding to the constraint (jlOp is sufficient 
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to produce the lowest perturbation we are interested in, yielding the perturbed mKdV 

Vt - Vxxx - Qv'^Vx = Wx{t, X), (19) 

Wxx-M(?{i)-wl)'^ (20) 

To see the effect of deformation more closely we construct explicitly 1- soliton solution, referring to [20] 
for the details of the the details on the general N-soliton. The accelerating soliton of the perturbed 
mKdV, shown in Fig. 2, has the explicit form 

v{x, t) = — sech^, ^ = k{x — vt) + <P, v = vq + Vd (21) 

where vq = Ak? is the usual constant velocity of the mKdV soliton, while Vd = is the unusual time- 
dependent part of the velocity, induced by the deformation. Note that the perturbing function w{x,t) 
itself takes a self-consistent solitonic form 



Wx = c(t)Ksech(^tanh^ 

and drives the field soliton (j21|) to an accelerated motion again through its asymptotic value w{x,t)\^ 
c{t), sitting as a forcing term at the space boundaries. 



(22) 




Figure 2: Exact accelerated soliton in the (x,t)-plane for the integrable perturbed mKdV equation. 
. C. Integrable perturbation of the NLS equation 

For constructing the integrable perturbation of the NLS equation the stage is already set in Q, since 
with a reduction r = q* and a minimum deformation (jlOp . it would yield the NLS case 



iqt - Qxx - 2|g| q = g, 



(23) 



with the perturbative function g{x, t) subjected to a nonholonomic differential constraint 



Qx = -2iaq, = i{qg* - q*g) 



(24) 
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where = 5) G^2i = 5*' ^ii — ~^22 — ^ deforming functions coupled through the basic 

fields q, q* . The set of constraints (I24p may be simplified to a single differential constraint 

Hg) = gxxQ - gxQx - '^q^iqg* - q*g) = o. (25) 

Eliminating the deforming function g from (123p and (j25p we can further derive a new 4th-order NLS 
equation expressed through the basic field as 

{qdL - qxd^){iqt - qxx - 2\qW) + 2q\i{\qW + {qql - q*q^),) = 0. (26) 

As already mentioned, this perturbed system can be solved exactly through the ISM. Remarkably how- 
ever, we can actually avoid the involved procedure of the ISM in all such integrable nonholonomic 
deformations and can obtain the explicit soliton solutions for these deformed equations by just suit- 
ably deforming the known undeformed solutions. As a rule the perturbing affect keeps the form of the 
soliton intact, while it deforms some parameters like velocity, frequency etc., making them in general 
time-dependent . 

More precisely, in case of the deformed NLS, the original soliton velocity vq and the frequency ujq of 
its enveloping wave are changed with an addition of a deformed velocity Vd{t) = and the deformed 
wave frequency iOd{t) = —2^^j^, where Ai = ^ -|- irj. The function c{t) responsible for such deformation 
is linked to the asymptotic value of the perturbing function. The accelerating NLS soliton qualitatively 
looks like that of the mKdV when the dynamics of \q\ is plotted (see Fig. 2). 

We can generate again a two-fold integrable hierarchy for the perturbed NLS system. The first one 
is the known hierarchy of higher order NLS equations with a perturbation constrained by the same 
nonholonomic deformation (j25p . while the second one is a new integrable hierarchy, where the same 
perturbed NLS equation (j23p is perturbed by a hierarchy of higher order deformations of the form (jl2p . 
For example, the next to the lowest order deformation (j25p would include another deforming function 
e{x,t), in addition to function g{x,t) . Interestingly in this double-deformation the main equation would 
remain same as the perturbed NLS equation (p3|) . while the constraint equation for g would be changed 
to 

L{g) + 2i{eqx - qe^) = 0, 
or gxxq - gxqx - 2q'^{qg* - q*g) +2i{eq^ - qe^) = 0, (27) 

coupling to both the field q and the function e. The second deforming function e in turn is constrained 
exactly by the same differential constraint psp as 

L(e) = 0, or e^xq - exqx - 2q'^{qe* - q*e) = 0. (28) 

Even though these new integrable perturbed NLS equations seem to be rather academic, it is intriguing 
that, a model like our perturbed system (j23ll25p is implemented already in doped fiber for efficient optical 
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communication and therefore the implementation of the higher order integrable deformation discovered 
here for a a similar multi-doped system seems to be a promising possibility |21j . 

D. Integrable perturbation of the SG equation 

Notice that in the line of the present approach, the well known SG equation in the light-cone co- 
ordinates, which can be obtained through reduction q = r = 6^ from the AKNS system, can itself be 
considered as a deformation 9xt = —ig of the linear wave equation O^t = 0, with a constraint (jlOp on the 
deforming function g. This constraint can be easily resolved as = g = ism26, yielding the standard 
SG equation: 9xt = sin 20 . The situation however becomes more interesting if in the same equation 
9xt = —ig, the deformation g = ig is coupled to another deformation G^^^ , yielding the nonholonomic 
constraint (llip in the form 

gx = -2a9x + 2e, = 29xg, Cx = 0, (29) 

with the matrix elements Gy^ = a, = —G21 = ig, G^2 — ^21 ~ ~ 0' where e(t) is an 

arbitrary function of t . Expressing (I29p through the perturbing function g we obtain the nonholonomic 
deformation of the SG equation as 

9xt = g, (30) 
gxx9x - gx9xx + ^g9l + 2e{t)9xx = (31) 

Note that solving the constraint ()29p or equivalently (|3ip one can get different deformations of the SG 
equation. We derive one such interesting solution as 

g = e{t){a sm.29 + 13 cos29), 

a = e{t) {a cos 29 - (3 sin 29). ax = sin 29, Px = cos 29, (32) 

which yields a new integrable deformation of the SG as 

9 xt = e{t) {a sin 29 + 13 cos 29), (33) 

with a, f3 as defined in (I32p . It is intriguing to note that, the new SG equation (j33p has an additional 
cos 29 part along with the traditional sin 20 term and the coefficient ea, which usually corresponds to the 
particle mass becomes a space-time dependent function. It is worth mentioning that, though at e = one 
can recover from (j29p the undeformed standard SG equation, as shown above, we can no longer obtain 
the standard SG equation from its deformation ()33p at e ^ 0. This new equation ()33p . which is related 
to the system considered in [16], should however be examined carefully for the applicability of the ISM. 
Nevertheless, assuming the existence of its kink solution we present in Fig. 3 the accelerated kink for the 
field 9 and the soliton solution for the perturbing function for the integrable perturbed SG equation. 

Continuing the process of higher order differential constraints, by introducing more and more terms 
with \~'^,n = 3,4,... in V^'^'^-^^X), we can generate a hierarchy of integrable higher nonholonomic 
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Figure 3: Exact soliton solutions for the integrable perturbation of the sine-Gordon equation, a) Acceler- 
ating kink solution for 9{x, t). b) Perturbation g{x, t) in the solitonic form with nontrivial time-dependent 
asymptotic c{t) = 0.5t 

deformations in most of the cases. In some cases, as the explicit calculations show, one has to add more 
than one higher terms to get the first nontrivial result. 

VI. PRESENT NONLINEARIZATION SCHEME & INTEGRABLE PERTURBATIONS 
IN THE BACKGROUND OF EARLIER RESULTS 

We compared our nonlinearization scheme for generating integrable systems with the AKNS construc- 
tion in Sec. III.B. It is worth mentioning here that the idea of this nonlinearization bears some analogy 
as well as differences with the well known Zakharav-Shabat (ZS) dressing method [22j. In the dressing 
method the Lax pair as well as the soliton solution are constructed starting from the asymptotic vacuum 
solution. In our nonlinearization scheme on the other hand the Lax operators are build up starting from 
the linear field equation, based only on the dimensional argument and the notion of the building blocks. 
In analogy with the ZS dressing method it is however tempting to construct in our scheme the soliton 
solution for the nonlinear equation starting from the solution of its linearized equation. 

Extending our construction further by nonholonomic deformation, we have gone beyond the standard 
AKNS method and considered perturbing matrices of scaling dimension L > in the Vn^X) operator. 
This amounts to extending this operator as Vn^mW = Vn{X) + V^^^'^^A) by considering, together 
with the positive powers of A upto N, all its negative powers upto M. We have shown through such 
deformations a beautiful universality of the recently found nonholonomic deformation of the KdV system 
OEj, by extending this concept to all members of the AKNS family, including the mKdV, the NLS and 
the SG equations. At the same time we have discovered a two-fold integrable hierarchy for them, with the 
explicit finding of some interesting multiple deformation. Our construction allows the application of the 
exact ISM for obtaining N-soliton solutions for both the field and the perturbing functions, which shows 
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an unusual accelerated motion for the solitons. However a significant advantage of such deformations is 
that, we can actually avoid the application of the ISM in these cases and construct the exact soliton of 
the deformed equations by deforming the well known soliton solutions. Such perturbations only deforms 
the time evolution of the spectral data, which deforms in turn the soliton velocity and for the complex 
fields also the frequency of the enveloping wave. 

It is true that in the long development of the theory of integrable systems, the use of negative powers 
of A in the Lax operator V{X) was implemented also in few earlier occasions. However a consistent, 
systematic and full utilization of this procedure including the simultaneous presence of higher positive 
and negative powers of A, associated to all equations in the AKNS family, and most importantly, the 
possibility of introducing new perturbing functions as multiple deformations, which are our main emphasis 
here, have never been undertaken. Similarly, the existence of a two-fold integrable hierarchy and exact 
N-soliton solutions for both the field and the perturbing functions in all integrable perturbed systems 
of the AKNS family, as well as the possibility of boosting the soliton to accelerated motion by the 
perturbing function through its time-dependent asymptotic at the space boundaries, as revealed here, 
were not explored or clarified in any earlier occasions. 

As we understand, the AKNS themselves have considered integrable SG and Maxwell-Bloch systems 
going upto A~"^ [Ij. A -ve hierarchy for the mKdV model was also considered in fl6j. However these work 
either did not consider both -|-ve and -ve powers of A simultaneously, which would result to the integrable 
perturbations, or they went away from the perturbative models by considering only the equations for the 
basic field. Thus these results did not come close to ours, though we can derive them as particular cases 
of our more general integrable model. 

However an important physically applicable model considered earlier [23j can be represented by the 
lowest order deformation of the NLS equation (|23ti25p found through our construction. This in turn 
indicates the physical significance of our integrable perturbed equations and opens up the possibility 
for application of the hierarchy of perturbed NLS equation, found here, to the nonlinear fiber optics 
communication through media with multiple doping, by switching to higher order deformations ^21j. 

We should pay special attention to another class of interesting source equations proposed by Melnikov 
|17j and show that our integrable perturbed models stand as complementary to the Melnikov's models 
and can make contact with them at a highly degenerate limit. In the Melnikov's formalism a set of 
eigenvalues A„, 1, 2, . . . , appear explicitly, which are needed to be distinct and strictly nonvanishing for 
any construction of those models. However as we show below, our construction could be related in fact 
to the complementary limit of the Melnikov's systems, when all these eigenvalues become degenerate and 
moreover go to zero, simultaneously. Miraculously, Melnikov's system, in general too complicated, sim- 
plifies drastically at this highly degenerate limit and reduces to our simple perturbative model, allowing 
exact accelerating solitons and suitable for physical applications. 
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Let us consider as a demonstration Melnikov's source equation for the NLS system given as 

iut + u^^ + 2\u\'^u = Y,i<l>l + i'l) (34) 

n 

and its complex conjugate, together with a eigenvalue problem at discrete eigenvalues An, n = 1, 2, . . . , A'' 
for a set of complex functions +ipn)- 

(pn,x + Ulf^n = >^n4>n 
1pn,x + U*(j)n = -Kipn, (35) 

and their complex conjugates. As such this set of 2(1 + 2A^) number of coupled complex equations is 
evidently much more complicated than our perturbed NLS (j23ti25|) and its application to physical models 
seems to be rather obscure. If however we denote the rhs of by J^ni'Pn'^'^n) = 9 ^-nd the combination 
Yl,n{4>ni^n — 4'n''Pn) — ^) then after some algebraic manipulation from ([55]) we find 

n 

b^ + ug*-u*g = (36) 

It is clear that the system closes and immensely simplifies only when all eigenvalues become the same 
as well as vanishing: A„ = 0, n = 1, 2, . . . , A^, which however is opposite to the Melnikov's theory, which 
strictly demands all 7^ 0. We note immediately that at this vanishing limit from (j36p we get 

5a; = -2ub, bx = u*g - ug*, (37) 

which generates our perturbed NLS (|23ti25p and for the real field reduction u* = u, g* = —g reproduce 
further the perturbed mKdV equation (j20p . Therefore we may conclude that at highly degenerate limit 
and in a situation complimentary to the Melnikov's assumption, the well known source equations re- 
duce to the lowest deformation cases of our perturbed equations. However the possibility of generating 
higher order deformation for perturbing functions constituting the novel integrable two-fold hierarchy of 
integrable perturbed equations, we have discovered, seems to be absent in the Melnikov's formalism. 

VII. CONCLUDING REMARKS 

The present scheme of nonlinearizing linear equations to integrable systems, though bears similarity 
with the AKNS approach, differs in its motivation, simplicity and generality. Our procedure does not start 
from a given integrable nonlinear equation, as customary in the AKNS method, but aims to construct it 
by nonlinearizing a given linear field equation and is applicable without much difficulty to the AKNS as 
well as to the KN family. Moreover due to identification of the building blocks of the Lax operators and 
the effective use of the physical notion of dimensional analysis, our construction of the crucial time-Lax 
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operator becomes much simpler. In place of 4(A^ + 1) unknown functions, in general complex, 

in the AKNS scheme obeying coupled partial differential equations, we have only + 1 unknown real 
integers to be determined from simple algebraic equalities. This simplicity of our scheme becomes more 
evident at higher values of A^. 

The extension of the notion of integrability from the KdV equation to the whole family of the AKNS 
system was a major breakthrough achieved in early seventies. We now find a similar extension for the 
nonholonomic deformation of the KdV equation, discovered recently [6, 7J, to all members of the AKNS 
family, which shows the universality of such integrable deformations beyond a single known example 
[6]. Moreover, by discovering an integrable hierarchy of increasingly higher order deformation for each of 
these cases, we establish that this universality, in fact is two-fold. One is the known AKNS hierarchy with 
an additional fixed deformation, while the other is a new integrable hierarchy for each of the NLS, the 
KdV, the mKdV and the SG equations, deformed by perturbing functions with increasingly higher order 
nonholonomic constraints. Interestingly, in the second scenario nonlinear integrable equations can be 
generated by perturbations only, even from a linear dispersionless or trivial equation like Qt = Qx, Qt = 0. 

Such integrable perturbations are achieved through an extension of our nonlinearization scheme, where 
we go beyond our building block conjecture and introduce new perturbative functions at every higher 
step of deformation of the time-Lax operator, making the fullest use of its negative power expansion in 
the spectral parameter. It seems that, in spite of numerous work in this well studied subject, these two 
simple steps were not considered earlier consistently, thereby missing a whole class of integrable perturbed 
equations, which we discover here. These integrable perturbed equations with nonholonomic constraints 
are different and much simpler than the traditional source equations of Melnikov. Interestingly however, 
they are related in a exceedingly subtle way, being complementary to each other and linked at a highly 
degenerate limit. Nevertheless such a relation can be shown to exist only at the lowest order deformation 
of our equations, while the novel hierarchy of integrable deformation, that we find for our perturbed 
systems, is totally absent in Melnikov's source equations. 

The nonholonomic deformation inducing new integrable equations presented here are exactly solvable 
through the inverse scattering method. However an important advantage of such deformed integrable 
equations is that, one can completely bypass the complicated ISM for extracting their exact soliton 
solutions and construct them by merely deforming the well known solutions of the unperturbed equations. 
One finds that under such perturbations the form of the solitons remain the same, while the soliton 
velocity and the wave frequency get additional terms, capable of introducing acceleration or deceleration 
in the solitonic motion. 

The explicit soliton solutions for the simplest integrable perturbed NLS, mKdV, KdV and SG equa- 
tions, presented here show such unusual accelerated motion forced by the perturbation. Within this 
framework of integrable perturbation, one can also arrange these systems to receive consistent multiple 
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forcing, which are expected to have significant physical appHcations, especially in fiber optics communi- 
cation in doped media [HI [231 121] • 

Extentions of the nonholonomic deformation to the KN family of equations [23] as well as to 2 + 1 
dimensional models are exciting future problems. Another challenging problem, which we could not solve 
here, is to construct the soliton solution of the nonlinear equation from the solution of its linearized 
equation, by using the present nonlinearization scheme, which might define some new dressing method. 
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